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NMR investigation of contextuality in a quantum harmonic oscillator 
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Physical potentials are routinely approximated to harmonic potentials so as to analytically solve 
the system dynamics. Often it is important to know when a quantum harmonic oscillator (QHO) 
behaves quantum mechanically and when classically. Recently Su et. al. [Phys. Rev. A 85, 
052126 (2012)] have theoretically shown that QHO exhibits quantum contextuality (QC) for a 
certain set of pseudospin observables. In this work, we encode the four eigenstates of a QHO onto 
four Zeeman product states of a pair of spin-1/2 nuclei. Using the techniques of NMR quantum 
information processing, we then demonstrate the violation of a state-dependent inequality arising 
from the noncontextual hidden variable model, under specific experimental arrangements. We also 
experimentally demonstrate the violation of a state-independent inequality by thermal equilibrium 
states of nuclear spins, thereby assessing their quantumness. 
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I. INTRODUCTION 

Quantum contextuality (QC) states that the outcome 
of the measurement depends not only on the system and 
the observable but also on the context of the measure¬ 
ment, i.e., on other compatible observables which are 
measured along with [1-4]. 

Let us consider a pair of space-like separated entangled 
particles, with local observables A and C belonging to the 
first particle, and B and D to the second. We assume 
that these observables are dichotomic (i.e., can take val¬ 
ues ±1) and that the pairs ( B,C ), ( C,D ), and 

( D,A ) commute. 

Classically, one assigns objective properties to the par¬ 
ticles such that D behaves identically on the state of the 
system irrespective of whether it is measured in the con¬ 
text of A or in the context of C, even though A and C 
are not compatible [5, 6]. Such measurements are said to 
be context independent. Classically, one can pre-assign 
values (a, c) to ( A , C) of the first particle independent 
of the measurement carried out on the second particle. 
Similarly, for the second particle one can pre-assign val¬ 
ues (&, d) to ( B , D) independent of the measurement car¬ 
ried out on the first particle. In these pre-assignments, 
implicit is the assumption of noncontextual hidden vari¬ 
ables, which predict definite measurement outcomes in¬ 
dependent of measuring arrangement. If we pre-assign 
values to observables such that A,B,C,D = ±1, it fol¬ 
lows that AB+BC+CD—AD = ±2 and the expectation 
value, 

I = (AB + BC + CD — AD) 

= {AB) + {BC) + {CD) - {AD) < 2 (1) 

[7]. This inequality often known as CHSH inequal¬ 
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ity arises from noncontextual hidden variable (NCHV) 
model and must be satisfied by all classical particles. 

Now let us see the implication of the quantum the¬ 
ory. Let Alice and Bob share a large number of singlet 
states: (|01) — 110 ))/a/ 2 = — (|H—) — |—|-))/v^2, where 
|0) and |1) are eigenkets of Pauli -0 operator (a z ) and 
|±) = (|0) ± \l))/y/2. Alice measures on her qubit ei¬ 
ther erf or erf, while Bob always measures erf. Let us 
compare the results of only those measurements in which 
Alice has obtained outcome +1. If Alice measures erf, 
then Bob’s qubit collapse to |1) = (|+) — |— ))/V2. In this 
context (i.e., erf), Bob will get both outcomes ±1 with 
equal probability. On the other hand, if Alice measures 
erf on her qubit, then Bob’s qubit collapses to |—) and in 
this context (i.e., erf), Bob will always get the outcome 
— 1. Hence the context dependency. 

Here we experimentally investigate QC of a quan¬ 
tum harmonic oscillator (QHO). There are a variety of 
quantum systems whose potentials are approximated by 
QHO. Consider for example the quantized electromag¬ 
netic field used to manipulate a qubit in cavity quantum 
electrodynamics [8]. Recently, QC in QHO has been the¬ 
oretically studied by Su et al [9] by mapping four low¬ 
ermost QHO states onto four pseudospin states. Such 
states can be encoded by qubit states, and QC can be 
studied by realizing the measurements of appropriate ob¬ 
servables. We realize this study using a nuclear magnetic 
resonance (NMR) quantum simulator [10] . 

In the following section we shall revisit the formula¬ 
tion of Su et al, and in section III, we describe the ex¬ 
perimental demonstration of state-dependent and state- 
independent QC using an NMR system. Finally we con¬ 
clude in section IV. 


II. THEORY 

Hong-Yi Su et.al. [9] have theoretically studied QC of 
eigenstates of 1D-QHO by introducing two sets of pseudo- 
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spin operators, 


r = (r X) r w ,r,) ) r' = (r;,r' r'j 


with components, 


T x = cr x <g> l,r y = <r~ (g> <Ty,r z = -a y <8> cr y , 

r^, = cr x g) a z , r' y = 1 ® (Jy, Y' z = —a x ® <r x . (2) 

where 1 is 2 x 2 identity matrix. Using these operators 
they defined the following observables, 

A = Y X = a x ® 1, 

B = T' X cos (3 + Y' z sin j3 = a x g (tr z cos ft — <j x sin /3), 

C r z (7 y g (7y 1 

D = Y' X cos 77 + Y' z sin 77 = a x g (er~ cos 77 — a x sin if). (3) 

The products which form the inequality expression (Eq. 
1 ) are 

AB = 1 g (cos /3 a z — sin /3 a x ), 

BC = —a z g (cos P u x + sin /3 a z ) 

CD = —a z g (cos 77 <j x + sin 77 a z ), 

DA = 1 g (cos 77 a, — sin 77 a x ). (4) 

Following commutation relations hold: [T,, F)] = 

0 (■ i,j = x,y,z), [Y x ,Y y ] = 2i Y z , [T(,, Y' y ] = 2 i Y' z and 
cyclic permutation of x,y,z. A,B,C,D have eigenval¬ 
ues ±1, with (A,B), (B,C), ( C,D ) and (D,A) forming 
compatible pairs. A, B , C and D are Hermitian, hence 
observables. One can verify that they are also unitary 
operators. Hong-Yi Su et.al. [9] have shown that, 


1 ( \i)qho = 2 ^ > 2 ’ when ’ (Z 3 ’ 7 ?)' = 


f (-T! - / 4 ! -3 tt/4)o 
I (3tt/4,tt/4)i 
1 (tt/4, 3tt/4) 2 
[ (- 37 r/ 4 ,- 7 r/ 4) 3 


where, ljj^ o is the expression on LHS of inequality 

1, l = 0,1,2 and 3, and, \0)qhoA^)qhoA^)qho and 
|3 )qho are first four energy eigenstates of 1D-QHO. Thus 
QHO violates the inequality (1) for certain observables 
and thereby exhibits QC. 

It is well known that only certain two-particle states 
violate the CHSH inequality (1). As shown in [11, 12] 
factorable states always satisfy inequality ( 1 ) for local 
observables (observables of the form Pg 1 or 1 gQ [13]). 
With maximally mixed state (1 g 1 /4) the inequality (1) 
is satisfied even with nonlocal observables (observables of 
the form PgQ [13]) in eq. (3), which is obvious from the 
fact that all the products in eq. (4) are traceless. How¬ 
ever, if the initial state is nonfactorable, we can always 
find observables such that inequality ( 1 ) is violated [ 11 ]. 
Although the pseudospin states {100), |01), 110), |11)}, are 
factorable, they still violate ( 1 ) since the observables in 
eq. (3) are nonlocal. Thus, we observe that even when 
a system is in a nonentangled state, measurements of 
nonlocal observables may lead to violation of noncontex- 



FIG. 1. Moussa Protocol for finding out the expectation value 
of the joint observable X\ X 2 X 3 i.e. {X 1 X 2 X 3 }. Here Xi s are 
mutually commuting unitary operators. 


tuality inequality [14]. 

State independent QC: There exist stronger in¬ 
equalities obtained from NCHV models which is violated 
by all states, including separable or maximally mixed 
states. If the initial state is maximally mixed, entangle¬ 
ment cannot be created by measuring whatever observ¬ 
able (local or nonlocal). This shows that entanglement 
is not necessary in general even in a bipartite system, to 
exhibit QC. Hence we conclude that, QC is more funda¬ 
mental or general than entanglement. Any system whose 
Hilbert space has dimension > 2 exhibits QC [4]. Even a 
single spin -1 particle (where entanglement has no mean¬ 
ing as far as spin degree of freedom is concerned) can 
exhibit QC [15]. 


III. EXPERIMENT 
A. State dependent contextuality 

For experimentally studying eq. (1), we need: (i) a 
physical representation of first four energy eigenstates 
{|0)q77O: |1 )qho, |2 )qho, |3)qt7o} of 1D-QHO , and (ii) 
a way to find out the expectation values for opera¬ 
tors AB, BC , CD, and DA. We encode the first 
four energy eigenstates of 1D-QHO onto the four en¬ 
ergy eigenstates (under secular approximation) of a pair 
of spin- 1/2 nuclei precessing in external static magnetic 
field: {|00), |01), 110), 111)}, which we call Zeeman prod¬ 
uct states. In fact any four arbitrarily choosen en¬ 
ergy eigenstates of 1D-QHO and also their superposition 
states exhibit QC [9]. The circuit shown in Fig. 1 is 
called Moussa protocol [16], and is used to extract the 
expectation value of observables in a joint measurement. 
This has subsequently been generalized by Joshi et al [17] 
to unitary operators. 

The three qubits for this experiment were provided by 
the three 19 F nuclear spins of trifluoroiodoethylene dis¬ 
solved in acetone-D 6 . Fig. 2(a) shows the structure of 
trifluoroiodoethylene along with the Hamiltonian param¬ 
eters in Fig. 2(b). The effective 19 F spin-spin (T?j) and 
spin-lattice (Ti) relaxation time constants were about 
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FIG. 2. (a) Molecular Structure, (b) chemical shifts (diagonal 
elements) and J-couplings(off-diagonal elements) in Hz of tri- 
fluoroiodoethylene, and (c) pulse sequence for pseudo-pure 
state preparation. The amplitude and phase of the shaded 
pulse is written over the top of respective pulse, and unshaded 
pulses are (n) x pulses. During the J evolution 1/(2J 23 ), (n) x 
pulses refocus J 12 and J 13 evolutions and retain only J 23 evo¬ 
lution. Similarly others. Last row represents Pulsed Field 
Gradients(PFG), which kills transverse magnetization or in 
other words they make coefficients of I x and I y zero [18] 


0.8 and 6.3 s respectively. The experiments were carried 
out at an ambient temperature of 290 K on a 500 MHz 
Bruker UltraShicld NMR spectrometer. 


AP|000) =-(/lz + l2z + l3z + 2/i z /2z 

+2/22/32 + 2/12/32 + 4/12/22/32)- ( 8 ) 

The first spin, Fi, is used as an ancilla qubit, and 
other spins, F 2 and F 3 , as the system qubits (see Fig. 
1). The initial Hadamard gate on the first spin prepares 
P|+oo) ■ To measure (AB)|qo) , we apply the corresponding 
controlled operations A and B as indicated in circuit 1. 
The transverse magnetization of the ancilla qubit will be 
proportional to the expectation value (AB)|qo)- The ab¬ 
solute value of (AB)|qo) is estimated by normalizing the 
value obtained in the above experiment with that ob¬ 
tained from a reference experiment having no controlled 
operations. Similarly we can measure the other expec¬ 
tation values ( 5 ( 7 ) 100 ), ((7-D)|oo), an d (AD)|oo), and de¬ 
termine the value of Iq. Other values I; are obtained 
by preparing the corresponding pseudopure states P|+oi), 
Pl+io), an< 4 pi+ 11 ) and applying the circuit 1 in each case. 

In our experiments, all the controlled operations were 
realized by numerically optimized radio frequency (RF) 
pulses obtained using GRAPE technique [20]. Each pair 
of controlled operations in circuit 1 was realized by a 
GRAPE sequence with a duration of about 23 ms (having 
RF segments of duration 5 ps) and an average Hilbert- 
Schmidt fidelity better than 0.99 over 10% variation in 
RF amplitude. 

We estimated the values for I; ( 1 ), for all the four 
eigenstates and independently varied both /3 and ?) over 
the range [— 7 r, 7 r] with increments of 7 r/ 4 . The re¬ 
sults are shown in Fig. 3. The maximum theoreti¬ 
cal violation is 2\/2 = 2.82. The experimental val¬ 
ues of maximum violations for Io, Ii, I 2 , and I 3 are 
2.40 ± 0.02, 2.45 ± 0.02, 2.39 ± 0.02, and 2.42 ± 0.03 
respectively. 


The thermal equilibrium state for the three spin sys¬ 
tem in the eigenbasis of total Hamiltonian (under secular 
approximation: { 1000 ), | 001 ),...}) is 

Peq = y + e (6) 

i -1 

where, Is is an 8 x 8 identity matrix, are spin an¬ 
gular momentum operators, and the purity factor e = 
fiiuo/(8kT) is the ratio of the Zeeman energy gap to the 
thermal energy [18]. Unitary operation has no effect on 
the identity part, but modifies only the traceless devia¬ 
tion part. By applying a series of unitary and nonunitary 
operators (pulse sequence shown in Fig. 2 [19]), it is pos¬ 
sible to transform the equilibrium state to a pseudopure 
state 

Ppps = (1 — 4- e|000)(000| = — + eAp| 000 ) (7) 

which is isomorphic to the pure state |000) [10]. In the 
pseudopure state, the traceless deviation part has the 


B. State independent contextuality 

Hong-Yi Su et.al. also studied the state independent 
contextuality [9, 21]. They considered the inequality 
(arising from NCHV model) 

(P11P12P13) + (P21P22P23) + (Pn P32P33) 

+{Pi\P2iP3i) + (P 12 P 22 P 32 ) — ( P 13 P 23 P 33 ) < 4 (9) 

where Pij are the elements of the matrix P, 

/ 17 r' r 2 r'\ 

p= r; v x r x r; . (io) 

V r *r' r*r' z r y r; j 

Operators in each row of the matrix P commute with 

each other. Similarly in each column. P\j are dichotomic 
observables with measurement outcomes ±1. We can ver¬ 
ify the inequality 9 by preassigning the values ±1 to each 
observable P,; ? . 
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FIG. 3. I () ,Ii,l 2 ,and I 3 represents evaluation of expression 1 for eigenstates |0}Qno, |1 )qho, |2 )qho, and |3}qho respectively. 
The curved surface represents theoretical values, and the points are experimental values. Flat planes at 2 and —2 represents 
classical bound. 


Now introducing the operators from expressions (2), 
we find that the product of each row of the matrix P is 
identity (i.e., P 3 \ PfiPj-j, = 1) having eigenvalue +1. Sim¬ 
ilarly, the products along each of the first two columns is 
again identity. However, the product along the last col¬ 
umn, i.e., P13P23P33 = — 1 having the eigenvalue —1. No 
preassignment of ±1 to the various elements of P can sat¬ 
isfy the condition that, product along each row and along 
the first two columns be +1 and along the last column be 

— 1. This shows that quantum theory is not compatible 
with NCHV model. Further, for the above choice of op¬ 
erators, the expectation values for the first five operators 
in expression 9 are all +1 while that of the last term is 

— 1. Therefore, quantum bound for lhs of expression 9 is 
6 independent of initial state of the system. 

To investigate state-independent QC, we need to mea¬ 
sure joint expectation values of three observables. We 
again use the circuit 1 for this purpose. Taking advan¬ 
tage of the state independent property of the above men¬ 
tioned inequality, we choose thermal equilibrium state 
(6) as initial state. A (tt/2) v pulse was applied on the 
first spin to prepare the ancilla in a superposition state. 
Then state (6) takes the form: (1 — 4e)l s /8 + e(|+)(+| G> 
1 ® 1) + c(^2z + l3z)- 

All the controlled Pij operations were realized using 
the GRAPE sequences having average fidelities better 
than 0.99 over 10% variation in RF amplitude. The total 
duration of the RF sequences (for each term in 9) were 
about 40 ms. Experimentally obtained value of lhs of 
inequality 9 is 4.81 ± 0.02. 

Thus we observed a clear violation of the classical 
bound. However it is still lower than the quantum limit. 


The reduced violation is attributed to decoherence (T 2 
decay) and imperfections in RF pulses. 

IV. CONCLUSION 

We have experimentally demonstrated the quantum 
contextuality exhibited by first four energy eigenstates of 
a one dimensional quantum harmonic oscillator through 
the violation of an inequality obtained from a noncontex- 
tual hidden variable model. The continuous observables 
of the harmonic oscillator are mapped onto pseudospin 
observables which are then experimentally realized on a 
pair of qubits. We used Moussa protocol to retrieve the 
joint expectation values of observables using an ancil¬ 
lary qubit. Our quantum register was based on three 
mutually interacting spin-1/2 nuclei controlled by NMR 
techniques. We also demonstrated a violation of an in¬ 
equality formulated to study state-independent contex¬ 
tuality, by measuring a set of expectation values on the 
thermal equilibrium states of the nuclear spins. The re¬ 
sults of the experiment not only establish the validity of 
quantum theoretical calculations, but also highlights the 
success of NMR systems as quantum simulators. 
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